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The dynamics of the first order spatial coherence, (? , for a polariton system in a high-Q GaAs 
microcavity was investigated on the basis of Young's double slit experiment under 3 ps pulse excita- 
tion at the conditions of polariton Bose-Einstein condensation. It was found that in the process of 
condensate formation, the coherence buildup time increases almost linearly with the distance indi- 
cating that the coherence expands with almost constant velocity of about 10 8 cm/s. The coherence 
is strongly affected by exciton reservoir and polariton decay kinetics. It is smaller than coherence 
in thermally equilibrium system during the growth of condensate density and well exceeds it at the 
end of condensate decay. 



One of the most important characteristics of Bosc- 
Einstein condensate is the spatial coherence or the off- 
diagonal long range order, i.e. the property of the system 
to share the same wave function at different points sep- 
arated by a distance larger than the thermal de Broglic 
wavelength. To understand the processes governing the 
Bosc-Einstcin condensation (BEC), it is important to 
know how fast the coherence is established throughout 
the system during the condensate formation. This ques- 
tion was addressed theoretically in Refs. 0, where it 
was shown that in the process of a condensation particles 
first relax to the low energy or so-called coherent region 
where the kinetic energy of the particle is of the order of 
its interaction energy with other particles. Second, the 
fluctuations of density are smoothed out and a "quasicon- 
densate" is formed. Further, the phase fluctuations dis- 
appear resulting in the long range order formation within 
the system signifying the onset of the "true condensate" . 
The timescales of these processes are quite different: the 
relaxation to the low energy region is mainly determined 
by the stimulated scattering processes, whereas the qua- 
sicondensate formation is governed by the intcrparticlc 
interactions. Experimentally, the dynamics of spatial co- 
herence formation was first studied for a gas of ultracold 
atoms in Ref. Q, where it was found that the coherence 
expands with a constant velocity of about 0.1 mm/s. 

In this Letter we discuss the expansion of the spa- 
tial coherence in a condensate of mixed exciton-photon 
states, polaritons, in a semiconductor microcavity (MC) 
with embedded quantum wells. The bosonic statistics of 
these particles and the extremely light effective mass m 
(~ 1CP 5 of the free electron mass m e ) allow one to ob- 
serve MC polariton BEC up to the room temperatures, 
which inspired a considerable attention to this system 
in the last decade. After the first experimental demon- 
stration of MC polariton BEC [H a number of related 
phenomena in the MC polariton system have been ob- 



served and discussed: superfluidity quantized vor- 
tices , spin-Meissner effect and Josephson effect [H 
(see [9j] for a review). Compared with atomic BEC, the 
MC polariton condensation is highly specific: the system 
of MC polaritons is strongly nonequilibrium due to the 
short lifetime of the quasiparticles and the presence of 
pumping [Toj ] . 

The first order spatial coherence function, is re- 
lated to the off-diagonal elements of the coordinate de- 
pendent polariton density matrix g(p, p') as 
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and can be probed in the MC polariton system by mea- 
suring the interference of th e lig ht emitted from differ- 
ent points on the sample 0, 11, EH- Dynamics of 
in a process of MC polariton condensate formation was 
studied experimentally using a Michelson interferometer 
in [l3| in a CdTe-based MC for the fixed separation be- 
tween condensate regions Ax = 8.5 ^m and very recently 
in [3] in a GaAs-based MC, where the temporal dynam- 
ics of was related with the transition from photon to 
polariton lasing. A systematic study of spatial-temporal 
coherence expansion in a process of condensate formation 
is still missing to the best of our knowledge. 

In the present work we study the spatial coherence 
dynamics of the polariton system in a high-Q GaAs- 
based MC in the geometry of Young's double-slit ex- 
periment ll| for different separations between the slits 
under pulsed excitation. This allows us to trace the de- 
pendence g^(Ax) at different times and extract the dy- 
namics of the coherence length r c . Maximum r c w 17 fj,m 
is reached approximately simultaneously with the max- 
imal ground state population, at the excitation density 
P not strongly exceeding threshold density Pthr and r c 
decreases with P at P > Pthr- We have found that 
the coherence buildup time grows almost linearly with 



2 



distance Ax, that is the coherence expands with almost 
constant velocity. We show that formation of the spa- 
tial coherence under pulsed excitation with P < 2P t h r is 
governed by the processes of polariton relaxation towards 
the ground state, whereas amplitude and phase fluctua- 
tions of the quasicondensate are not suppressed in MC 
polariton system. 

The sample is a half wavelength MC with Bragg re- 
flectors made of 32 (for the top mirror) and 36 (for the 
bottom mirror) AlAs and Alg.13Gao.87As pairs. It has a 
Q-factor of about 7000 and the Rabi splitting of 5 meV. 
All the experiments were performed at 10 K and photon- 
exciton detuning of —9 meV. The sample was excited 
by the linearly polarized radiation of a mode-locked Ti- 
sapphire laser generating a periodic train of 2.5-ps-long 
pulses at a repetition rate of 76 MHz. The angle of in- 
cidence of the excitation laser beam was 60° with re- 
spect to the sample normal and the excitation photon 
energy corresponded to the reflection minimum of the 
mirror at this angle and was 11 meV above the energy 
of the bare exciton. The beam was focused in a 20 /im 
spot on the sample surface. The spot was imaged with 
magnification of T = 6 on the light-absorbing plate with 
two transparent parallel slits. The interference pattern 
of the emission coming from the regions of the sample 
selected by the two slits was formed on the Hamamatsu 
streak camera slit, like in Ref. ll|. Different slits pairs 
with separations X between slits in each pair ranging 
from 20 to 160 /im and slit width 5 /im was used to 
study g^(Ax), where Ax = X/T. Spatial coherence 
was extracted as the visibility of the interference pattern 

9 ^ (Imax Imin) / (^max ~t~ -^min) ; where Ifnin <md Imax 

arc minimal and maximal intensities within one period of 
interference pattern, averaged over all the observed peri- 
ods. The time resolution in this case was 3 ps. To study 
the energy-time dynamics of the MC emission it was dis- 
persed by a spectrometer coupled to the streak camera. 
In this configuration spectral and temporal resolutions 
were 0.25 meV and 20 ps, respectively. 

To convert the intensity I(t), measured by the streak- 
camera to the number of polaritons TV at states with 
wavevectors |fc| < 3 /im -1 , corresponding to the collec- 
tion aperture (above BEC threshold they mainly con- 
centrated at k ps 0), the integrated intensity of the MC 
photoluminescencc (PL) Ipl was measured by the sen- 
sitive power meter. Then the number of emitted pho- 
tons per laser pulse N pu i se was evaluated by the relation 
Np U i S e = 2Ipl/ (ffiuj), where factor 2 takes into account 
two directions of photon emission, / is the laser pulse 
repetition rate and hw is the energy of emitted pho- 
tons. Finally, the number of polaritons was evaluated 
by the relation TV(t) = N pu i se TLpI(t)/ J I(t)dt, where 
Tip = tc/\C\ 2 « 3 ps is the polariton lifetime at the 
bottom of the lower polariton (LP) branch (C is the pho- 
ton Hopfield coefficient, tc is the lifetime of photon in a 
cavity, which was determined from the linewidth 7c of 
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FIG. 1: Dynamics of the energy (full squares, left axis) and 
FWHM (empty circles, right axis) of the LP emission line. 
Emc and Elp are marked by black arrows. Inset shows the 
number of LPs near k — 0. Excitation power is l.8Pthr- 



strongly photonlike polariton: tc — h/jc)- 

At low excitation density P the PL dynamics of the LP 
branch is relatively slow, and the angular distribution of 
intensity indicates a bottleneck effect. As P is increased 
above the BEC threshold P t hr = 0.7 kW/cm 2 (this value 
corresponds to the time-averaged power of pulsed excita- 
tion), a fast and intense component in the PL dynamics 
corresponding to k ps appears. Onset of the fast compo- 
nent is accompanied by the blueshift of the spectral line 
and decrease of its width, which becomes determined by 
the spectral resolution (Fig. [T]) . The energy position of 
the spectral line is close the bare MC mode Emc just 
after the excitation pulse and relaxes to the energy of LP 
mode at low particle density Elp with time. The max- 
imum population (inset in Fig. [TJ is reached when the 
spectral line energy is between the MC and LP modes 
indicating that BEC is observed in the strong coupling 
regime despite the high particle density (IB- 17|. 

Below P t hr interference fringes in the double slit ex- 
periment are observed only for the smallest studied slits 
separation, and g^'(Ax = 3 /im) at 0.9Pt/ lr is less than 
0.3. Above Pthr the interference fringes are well resolved 
up to Ax = 20 /im. Figure HJa) shows the time dynam- 
ics of the interference of the emission coming from the 
sample regions separated by Ax = 3 and 17 /im. 

The dynamics of the spatial coherence g^ 1 for different 
separations Ax at P = 1.8P t hr is presented in Fig. [^b) 
together with the dynamics of the polariton number TV 
at the bottom of the LP branch. Figure [2jb) shows that 
the maximum value of g^ 1 ' decreases with increased Aa; 
and the decay of g^ in the whole range of Aa; occurs 
much slower than that of the condensate density. 

Interestingly, Fig. [UJb) shows that the coherence 
buildup time increases with an increase of Ax, indicating 
a finite velocity of coherence expansion. We define the 
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FIG. 2: (a) Streak camera images of the interference patterns of the emission coming from regions of the sample separated by 
Ax = 3 and 17 fxm. Horizontal axis corresponds to the angle of the emission passed through the slits with respect to the sample 
normal, (b) Dynamics of the spatial coherence for different Aa; (left axis) and number of particles near the bottom of the LP 
branch (right axis). Arrows mark the times to, 5, when g' 1 ' reaches half of its maximal value, (c) Dependence of the coherence 
buildup time to. 5 on Ax. Solid curve shows linear fit. The slope of the dependence yields the coherence expansion velocity u c . 
(d) Dependences of spatial coherence function on Ax at different times t. Experimental results are shown by symbols, solid 
lines show dependences calculated at the same times as experimental ones for the condensate size L — 17 /im and dashed line 
is dependence calculated at t = 160 ps and L = 10 ^m. Excitation power is 1.8Pthr- 



coherence buildup time to.s(Ax) as the time when g^ 1 ' 
reaches half of its maximum value for a given Ax. These 
times arc marked by arrows in Fig. [^b). Figure J^c) 
shows that the dependence of io.5 on Ax is close to lin- 
ear indicating that the coherence expands with almost 
constant velocity v c = 0.6 • 10 8 cm/s. 

Figure [3] shows the coherence length r c , which is de- 
fined from g^'(r c ) = 1/e and number of polaritons at 
the bottom of the LP branch N for P > Pthr- At the be- 
ginning of the BEC, r c grows almost linearly with time, 
reaches its maximum and decays afterwards. It follows 
from Fig. [3j that r c and N reach their maximal values 
almost at the same time, but the decay of N at high 
P occurs much faster than that of r c : after ~ 30 ps 
of relatively fast decay, in the dynamics of r c appears a 
slow component preserving the value of r c at the level 



0.5. . . 0.8 from its maximal value for several tens of pi- 
coseconds in the range where N decreases by more than 
an order of magnitude. Furthermore, Fig. [3] shows that 
the buildup of the coherence at P = 1.2P t / lr begins when 
the particle number N is more than one order of magni- 
tude smaller than that at P = 4.2P t /„. and that the max- 
imal value of r c decreases with P at P > 2P t h r in spite of 
a strong (about an two orders of magnitude) increase of 
the condensate density. These facts show that the length 
of coherence, as well as the rate of its increase, is not 
solely defined by the condensate occupation because of 
the nonequilibrium nature of the polariton BEC. 

In a 2D system of a finite size L under the conditions 
of thermal equilibrium the number of particles in the 
ground state can be estimated as 

N = N- N', (2) 
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FIG. 3: Dynamics of the coherence length (full squares, left 
axes, linear scales) and number of particles near the bottom of 
the LP branch (empty circles, right axes, logarithmic scales) 
at different excitation powers. 



where TV is the total number of polaritons (neglecting 
those in the reservoir) , TV' is the number of polaritons in 
all states but the ground, 



TV' 



kdkL 2 



/ L 

■(e- 
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(3) 
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where T is the temperature measured in the units of en- 
ergy. It is assumed in Eq. ([3]), that the chemical potential 
[i = 0. It follows from Figs.[2jb) and [3] that the onset of 
spatial coherence at P < 2P t hr starts at TV ~ 10 2 . This 
value is less than the estimated from Eq. © TV ss 400 at 
T = 10 K in the investigated LP system with the lateral 
size L rj 17 /im and m = 5-10 _5 ™ e . Thus, the formation 
of spatial coherence starts at negative chemical potential 
fi and hence governed by the relaxation process of polari- 
tons to the low energy region, while interaction-induced 
suppression of amplitude and phase fluctuations plays no 
role We note also, that at P < 2P t h r the estimated 

value of the interaction energy for the condensed polari- 
tons AE = aN Q /L 2 < 2 /teV (where a = 10~ 12 meVcm 2 
is the polariton-polariton interaction constant) is much 
less than thermal energy T « 1 meV. Hence, should 
be close to that for the classical noninteracting gas with 
the particle distribution function Nk [H| : 



gW(Ax) = 



% J™J (Axk)N k kdk + N 



L- 



f™ L N k kdk + N 



(4) 



where Jq(x) is the zero order Besscl function. 

The detailed modelin g o f the polariton distribution 
similar to that in Refs. [18J is beyond the scope of the 
present paper. Here we compare the experimental de- 
pendences g( l '(Ax) with ones calculated for the ther- 
mal Bosc distribution of the polaritons with the use of 
T = 10 K and the value of and chemical potential deter- 
mined from the measured polariton number iV in the low 
energy region. In the range of growing TV this simulation 
gives an upper limit for g^\ as the density of polaritons 
in excited states should exceed the equilibrium one, in 
order to provide the income of polaritons to the ground 
state. 



The results of calculation for three different stages: 
condensation onset, maximum condensate density and its 
decay, are shown in Fig. [2](d) for P = l.SPthr an d the con- 
densate size L = 17 /im determined from the experiment. 
It is seen that the experimental dependence g^(Ax) is 
slightly below the calculated thermally equilibrium one 
in the first stage (t = 45 ps), the difference strongly in- 
creases in the second stage (t = 85 ps) whereas in the 
third stage (t = 160 ps) experimental values g^(Ax) 
turn out well above the calculated ones. The letter is es- 
pecially surprising as it is difficult to expect the effective 
polariton temperature can be smaller than the bath tem- 
perature in the range of condensate decay. Measurements 
of polariton population along the LP branch (not shown) 
indicate that the distribution function approaches ther- 
mal one at the onset of condensation, but the occupation 
numbers of low k =/= states are slightly increased com- 
pared to thermal values, which explains a small discrep- 
ancy between the experimental and calculated g^\Ax) 
at t = 45 ps when TV w 350 < TV'. 

An increase in the discrepancy between the experi- 
mental and calculated g^(Ax) in the range of maxi- 
mum TV at t = 85 ps when TV w 7000 > TV' indi- 
cates that the polariton system becomes more noncqui- 
librium at high condensate density. The most probable 
reason for that is the runaway of condensed polaritons 
from the small bounded (~ 20 /im) photoexcited region 
due to their repulsive interaction with a dense exciton 
reservoir (ljj . In addition, for equilibrium system with 
TV 3> TV', the coherence is defined by the amplitude 
and phase fluctuations of the ground state wavefunction 
not taken into account in Eq. (j4|). Indeed, the quasicon- 
densate amplitude fluctuations are suppressed during the 
time ta ~ hL 2 /(aN) [l], determined by the interpar- 
ticle interaction. Taking a = 10~ 12 meVcm 2 , TV = 7000, 
L = 17 /tm we obtain value ta ~ 0.3 ns that exceeds the 
lifetime of polariton condensate at P = l.8Pth r , which is 
about 100 ps (Fig. [3]). While amplitude fluctuations and 
especially phase fluctuations are not suppressed in our 
experiment, the coherence length is determined by the 
healing length f = , hL , r 111, H , which is about 20 um 
at P = 1.8P t hr (N ~ 7000) and decreases up to 10 /jm 
at P = 4.2P t7ir due to increased TV (Fig.[5J. 

Finally, let us discuss the reason for the unexpect- 
edly high coherence in the decaying polariton system at 
t ~ 160 ps with TV w 200 < TV' which exceeds markedly 
the coherence in the thermally equilibrium system. Note 
that the calculations underestimate experimental values 
of g^ even if one takes unrealistically small conden- 
sate size L = 10 /tm (dashed line in Fig. [5U). At large 
t > 100 ps the reservoir becomes highly depleted, and 
the equilibrium in the system is established mainly via 
rather ineffective excitation of condensed polaritons by 
acoustic phonons with characteristic time exceeding the 
intensity decay time r. In this case the wave function at 
a given distance p decays as tjj(p) oc exp (— j-). It fol- 
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lows then, that both diagonal and off-diagonal elements 
of the polariton density matrix decay approximately with 
the same rate 1/t resulting in a nearly constant ratio of 
the polariton numbers in the ground state to that in the 
excited states well exceeding the decaying equilibrium ra- 
tio at long times. As a result weakly decreases with 
time even in the range where N < N', see Eq. (p}, in 
agreement with the experiment. 

To conclude, we have studied the dynamics of the spa- 
tial coherence for a condensate of MC polaritons under 
pulsed ps-long excitation for different excitation powers. 
It has been found that in the process of condensate for- 
mation, first order coherence expands with almost con- 
stant velocity of about 10 8 cm/s. We have shown that 
the coherence is determined not only by the density of 
low-fc polaritons, but also by the reservoir population. 
The coherence length r c and the number of low-fc po- 
laritons N reach their maximal values almost simulta- 
neously, afterwards r c decays much slower than TV. We 
have demonstrated that the onset of spatial coherence is 
determined by the narrowing of polariton distribution in 
fc-spacc rather than formation of the condensate phase, 
and at high excitation density coherence is limited by 
condensate amplitude fluctuations. The true condensate 
i.e. macroscopic occupation of the ground state with sup- 
pressed phase fluctuation is not achieved under ps-long 
pulsed pumping. 
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